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Abstract: We consider supersymmetric (SUSY) models wherein the strong CP problem is 
solved by the Peccei-Quinn (PQ) mechanism with a concommitant axion/axino supermul- 
tiplet. We examine i?-parity conserving models where the neutralino is the lightest SUSY 
particle, so that a mixture of neutralinos and axions serve as cold dark matter (aZ\ CDM). 
The mixed aZ\ CDM scenario can match the measured dark matter abundance for SUSY 
models which typically give too low a value of the usual thermal neutralino abundance, 
such as models with wino-like or higgsino-like dark matter. The usual thermal neutralino 
abundance can be greatly enhanced by the decay of thermally-produced axinos (a) to neu- 
tralinos, followed by neutralino re- annihilation at temperatures much lower than freeze-out. 
In this case, the relic density is usually neutralino dominated, and goes as ~ (f a /N)/m d . 
If axino decay occurs before neutralino freeze-out, then instead the neutralino abundance 
can be augmented by relic axions to match the measured abundance. Entropy production 
from late-time axino decays can diminish the axion abundance, but ultimately not the 
neutralino abundance. In aZ\ CDM models, it may be possible to detect both a WIMP 
and an axion as dark matter relics. We also discuss possible modifications of our results 
due to production and decay of saxions. In the appendices, we present expressions for the 
Hubble expansion rate and the axion and neutralino relic densities in radiation, matter 
and decaying-particle dominated universes. 
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1. Introduction 



The Standard Model of particle physics is beset by two major fine-tuning problems. The 
first occurs in the Higgs sector of the electroweak model, where quantum corrections to the 
Higgs mass push up to the energy scale associated with the cut-off of the theory, where 
new physics is expected to enter. The electroweak fine-tuning problem is elegantly solved 
by the introduction of weak scale supersymmetry (SUSY)]!]]. A consequence of weak scale 
SUSY is that supersymmetric matter states should exist at or around the electroweak scale, 
and ought to be detectable at the CERN LHC0. 

The second fine-tuning problem arises in the QCD sector. Here, t'Hooft's solution 
to the U(1)a problem Q via instantons and the 9 vacuum requires the existence of a CP 
violating term 

in the QCD Lagrangian[Q]. A second contribution to Lq arises from the electroweak sector, 
and is proportional to arg(detAi), where Ai is the quark mass matrixQ. Measurements 
of the neutron EDM tell us that the combination 9 + arg{detM.) = 9 must be < 10 _11 |(|. 
Explaining why 9 is so small is known as the strong CP problem. The strong CP problem 
is elegantly solved by the Peccei-Quinn mechanism |7|] and its concommitant axion a(8|. 
Models of an "invisible axion" with PQ breaking scale f a /N > 10 9 GeV (with N being the 
color anomaly factor, which is 1 for KSZVQ models and 6 for DFSZ[|lO| models) allow for 
a solution to the strong CP problem while eluding astrophysical constraints arising due to 



energy loss from stars in the form of axion radiation [11]. 



Of course, the SUSY solution to the electroweak fine-tuning, and the PQ solution to the 



strong fine-tuning are not mutually exclusive. In fact, each complements the other[12], and 
both are expected to arise rather naturally from superstring models (13]. In models which 
invoke both imparity conserving SUSY and the PQ solution to the strong CP problem, 
the dark matter of the universe is expected to consist of a mixture of both the axion and 
the lightest-SUSY-particle (LSP). Many previous studies have focused on the possibility 



of an axino a as the LSP|L4|, [15], 16], giving rise to mixed axion/axino cold dark matter 
(CDM) []l7|] . In this paper, we explore instead the possibility that the lightest neutralino 
Z\ is the LSP, thus giving rise to mixed axion/neutralino {aZ\) CDM. In the case of mixed 
aZ\ CDM, it may be possible to detect relic axions as well as relic neutralinos (as WIMPs). 

The case of neutralino CDM in the PQ-augmented MSSM has been considered previ- 
ously by Choi et al. [18]. In Ref. |fDj| ], the authors considered the case of neutralino dark 
matter where m« > ■ They presented approximate expressions to estimate the relic 
density of neutralinos fi^/t 2 . They found that neutralinos can be produced thermally as 
usual, but also that their abundance can be augmented by thermal production of axinos in 
the early universe, followed by axino cascade decays into the stable Z\ state. However, the 
neutralino abundance could also be diminished by two effects. The first is that even after 
neutralino freeze-out, the additional late-time injection of neutralinos into the cosmic soup 
from axino decay can cause a re-annihilation effect. The second diminution effect occurs 
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when late decaying axinos inject entropy into the early universe after neutralino freeze-out, 
thus possibly diluting the frozen-out neutralino abundance. 

In this paper, our goal is to present explicit numerical calculations of the relic abun- 
dance of mixed aZ\ CDM in SUSY models. To this end, we include several new effects. 

• First, we note that in the PQMSSM with a Z\ as LSP, the dark matter will consist 
of an axion/neutralino admixture, so we always account for the axion contribution 
to the total DM abundance. 

• Second, we account for the measured abundance of CDM as is given by the recent 
WMAP7 analysis |L9|: 

Q DM h 2 = 0.1123 ± 0.0035 at 68% CL. (1.2) 

We seek to establish under what conditions of model parameters the theoretical 
prediction for the relic abundance of mixed aZ\ CDM can be in accord with the 
measured value. 

• Third, we seek to establish whether, when fulfilling the WMAP measured abundance, 
the mixed aZ\ DM is dominantly axion or dominantly neutralino, or a comparable 
mixture. Such an evaluation is important for determining the relative prospects of 
axion and WIMP direct detection experiments. 

The remainder of this paper is organized as follows. In Sec. ^, we introduce two 
SUSY benchmark models- the first from the minimal supergravity (mSUGRA) model[20] 



in the hyperbolic branch/focus point region pl| (where neutralinos are expected to be of the 



mixed higgsino variety) and the second arising from the gaugino AMSB[[22j] model[23, 24] 
(where the neutralino is of the nearly pure wino variety). These benchmarks will be used 
for illustrative calculations of mixed aZ\ CDM. In Sec. ||, we calculate the axino decay rate 
into all possible modes, while in Sec. |I| we discuss the thermal production of axinos in the 
early universe. In Sec. || we compute the temperature at which axinos decay, and the 
temperature T e at which they might dominate the energy density of the universe. In Sec. || 
and in Sec. [?], we discuss the computation of the axion and the neutralino relic density. In 
the case of neutralino production, we pay special attention to the neutralino re-annihilation 
which takes place due to injection of neutralinos into the early universe from axino decay. 
In Sec. |U we present our algorithm for estimating the mixed aZ\ CDM abundance in 
SUSY models. In Sec. |9|, we present our numerical results, and in Sec. |10| we discuss the 
effect of including saxion production and decay in the calculation. In Sec. [ll], we present 
our conclusions. In the appendices, we present explicit formulae for the Hubble expansion 
rate, the axion abundance, and the thermal neutralino abundance for a universe that is 
either radiation-dominated (RD), matter-dominated (MD), or decaying-particle-dominated 
(DD). 

2. Two benchmark scenarios 

Supersymmetric models with mixed aZ\ CDM will always have a neutralino dark matter 
component, which is mainly enhanced by axino production and decay, but might also be 
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diminished by further annihilations and by entropy injection. In addition, there will always 
be an axion component to the CDM, which arises as usual via vacuum misalignment. In 
SUGRA-based models with gaugino mass unification, the lightest neutralino Z\ is usually 
bino-like, and annihilation reactions are suppressed by their p-wave annihilation cross sec- 
tions, leading to an overabundance of neutralinos |25|]. Since SUSY models with mixed aZ\ 
CDM will have to crowd both Z\s and axions into the overall relic density, we find that the 
most promising models to realize mixed aZ\ CDM are those leading to an uneterabundance 
of the usual calculation of thermally-produced neutralinos, such as models with either a 



higgsino-like or wino-like Z\ state[26]: such neutralinos annihilate via s-wave rather than 
p-wave reactions, and naively suffer a dearth of relic abundance. 

We will thus choose two benchmark cases to examine mixed aZ\ CDM. The first 



occurs in the hyperbolic branch/focus point region|21] of mSUGRA[20| and has a mixed 
higgsino-wino-bino Z\ state, while the second is taken from the gaugino-anomaly-mediated- 
SUSY-breaking model (inoAMSB)|24j], which leads to a nearly pure wino-like Z\ state. 



2.1 Mixed higgsino DM in HB/FP region of mSUGRA model 

The spectra for our first benchmark model, BM1, is generated within the minimal su- 
pergravity, or mSUGRA model, using the lsasugra[0 spectrum generator from Isajet 



v7.81[28]. The input parameters are taken as 



• (m , m 1/2 , A , tan/3, sign(/j,) = (4525 GeV, 275 GeV, 0, 10, +) 

with mt = 173.3 GeV. The Z\ has mass 87.9 GeV, while the calculated thermal abundance 
of neutralinos from IsaReD p9[| is Vt s ~ d h? = 0.05. Many sparticle masses and low energy 
observables are listed in Table [j]. Since the weak scale value of the superpotential /i term 
is only 137.2 GeV, the Z\ is of mixed higgsino-bino-wino type. 



2.2 Wino-like DM from gaugino AMSB 

The spectra for benchmark point BM2 of Table [l] is generated within the gaugino AMSB 
model which is expected to arise in string theories with moduli-stabilization via fluxes and 
a large volume compactification[^]. In the inoAMSB model, the parameters are taken as 

• m 3/2) tan/3, sign(p) = (50 TeV, 10, +) 

while too = Aq = at the GUT scale. The GUT scale gaugino masses take the AMSB 

form, and are given at m GUT as M% = 1065.3 GeV, M 2 = 161.4 GeV and M 3 = -484 GeV. 

This choice leads to a neutralino state Z\ which is nearly pure wino, with a thermal relic 

abundance of Q~ d h 2 = 0.0016 as indicated in Table 111. 

Zi u 

3. Axino decays in the mixed axion- neutralino CDM scenario 

Since we assume m„ > m^, an important element in the thermal history of the universe 
with mixed aZ\ CDM will be knowledge of the axino decay width. 
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Table 1: Masses and parameters in GeV units for the HB/FP and inoAMSB benchmark points, 
computed with Isajet7.81 using m t = 173.3 GeV. 



The axinogluon-gluino coupling is given by 



C ~ a ~" = 1 16tt(/1/JV) ^^^SaFavv (3.1) 



Evaluation of the a —> gg decay width gives a result in accord with Ref. [18], and is given 
by 

Also, the axino-bino-5^ coupling is given by 

^ B = i^^h^]6B^ (3.3) 

where ay = 5 ,/2 /47r, = d lL B v — d u B^ and C a yy is a model-dependent factor equal to 
8/3 in the DFSZ model, and 0, 2/3, 8/3 for heavy quark charges &q = 0, —1/3, +2/3 in 
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the KSVZ model. For numerical purposes, we will take C a YY = 8/3 since that case occurs 
in both models. We then find: 



n 2 fi2 2 a (»)2 / m~ \ 

r( ^ Zi+7H ^^ ml ,, 4) 



,2 



and 



2 

ffl; 



• < 



2 




(3.6) 



where v± is the bino fraction of Z$ in the notation of Ref. 0]. 

The axino partial and total widths are exhibited in Fig. [I] for a), the HB/FP model 
and b). the inoAMSB model. We see that at very low values of ~ Trig , only the decay 
a — > Zij is open, and r„ is very small. As mj increases, additional decay modes become 
allowed, and contribute to Tq. In frame b)., enjoys a large increase when the decay to 
Z27 opens up, since in the inoAMSB case, the Z2 is mainly bino-like. Once the decay to 
gg opens up, this mode is dominant. Decays to ZiZ are always subdominant to Z^, owing 
to the fact that the photon has a larger component than does the Z. 

4. Thermal production of axinos in the early universe 



If the reheat temperature after inflation, Tr, exceeds the axino decoupling temperature) 14] 

2 /n 1 \ 3 



then reheat occurred before decoupling which allowed the axinos to reach thermal equi- 
librium. Their number density at the time of decoupling is given in terms of the yield 
variable, Y = n/s, as 

Y~ = ^\ T - 135C(3) 1 (42) 

with C(3) ~ 1.202, g*(Tdc P ) is the effective number of relativistic degrees of freedom at 
temperature Td cp , and s = ^-g*T 3 is the entropy density of radiation 1 . 

In the other case, where Tr < Td cp , the axinos were never in thermal equilibrium in 
the early universe. However, they could still be produced thermally via radiation off of 



other particles in thermal equilibrium |l||, 3C]. Here, we adopt a recent calculation of the 



thermally produced axino yield from Strumia [31|: 



Yr = ^^10-"giF(g s )-^—(^-) , (L.:i, 



10 4 GeV \f a /N 



1 For simplicity we assume g,s(T) = g,(T), where g*s is the total number of relativistic degrees of 
freedom at T used to compute s. 
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inoAMSB: m y2 = 50 TeV, tan|3 = 10, A = 0, sgn{|i| > 1 




m- [GeV] 



Figure 1: Partial and total decay width of axinos versus 777,5 for BM1 in the mSUGRA model with 
(m ,m 1 / 2l A 0l ta,n /3, sign(mu)) = (4525,275,0,10,+) and for BM2 in the inoAMSB model with 
m 3/2 = 50 TeV, tan = 10 and fi > 0. We take f a /N = 10 12 GeV. 



with F(g s ) ~ 20<7^ l n ^") an d gs is the strong coupling constant evaluated at Q = Tr. 

If we suppose that each produced axino will cascade decay into the Z\ state, we can 
naively estimate the decay-induced neutralino abundance by 

777 ~ 

n| h 2 = -^nT p h 2 . (4.4) 

Zl m- a 

This estimate of the neutralino abundance from the decay of thermally-produced axinos is 
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Figure 2: Plot of Q| h 2 versus f a /N for various values of T R and = 50, 200 and 400 GeV 
(lower-to- upper) . 

shown in Fig. § as a function of fa/N for Tr = 10 4 , 10 7 and 10 10 GeV (green, blue and 
red curves, respectively), and for = 50, 200 and 400 GeV (lower to upper curves). We 
see that the neutralino relic density from axino decay can be enormous, and it typically 
dominates over the thermal neutralino abundance as calculated, for instance, from the 
mSUGRA model. From our naive estimate, we see that a tremendous overproduction of 
neutralinos is obtained from thermal axino production and decay except when considering 
parameter regions of high f a /N or low Tr where the axino production rate is suppressed. 

5. Axino domination of the Universe 

Once the decay width of the axino has been calculated as in Sec. ||, we may then calculate 
the axino lifetime or, alternatively, the temperature of radiation at the time scale when 
nearly all axinos have decayed, Tp>. This is achieved by equating the Hubble and the decay 
rates 

H(T D )=Ta, (5.1) 

which implies 



T D = v / r (i Mp/(7r 2 ^(T D )/90) 1 / 4 , (5.2) 

where Mp is the reduced Planck mass M P ~ 2.4 x 10 18 GeV. The temperature T D also 
corresponds to the temperature of radiation when entropy injection from axino decays is 
nearly finished. We assume an exponential decay law for a which implies that the axino 



decays only cause a slow-down in the rate of cooling of the expanding universe [32]. In this 
case, Tp) should not be interpreted as a "second reheat" temperature. 

As the universe- filled with axinos and radiation- expands and cools, T drops below 
and the axinos become non-relativistic. At that point, the axino energy density decreases 
as T 3 , while the radiation energy density continues to decrease faster as T 4 . At some 
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temperature T e , the energy density of axinos will dominate the universe if the axinos have 
not yet decayed, i.e. provided T e > Tjj. By equating the energy density of radiation 
Pr = ir 2 g*T 4 /30 with the energy density of axinos p a \T<ma = msXaS, we can determine T e 
to bepf 

4 

T e = -m a Y a . (5.3) 

Furthermore, we can estimate the ratio of entropy after to the entropy before axino decay, 
r, again where T e > Td, byp2|] 

S f 4m a Y a 

r = -f- ~ = T e /T D . (5.4) 

In Fig. |3|a)., we plot the values of Td, T e and the temperature of neutralino freeze-out 
Tf r versus m a for the benchmark BM1: the HB/FP region of the mSUGRA model. We 
take f a /N = 10 12 GeV, and show T e for T R = 10 6 and 10 10 GeV. For values of T D < 2 
MeV (gray shaded region), the parameters would be likely excluded because axinos would 
dump entropy after BBN has started. For Td > Tt r — the high m„ region- axinos decay to 
neutralinos before freeze-out. In this case, the neutralinos from axino decays thermalize and 
the neutralino relic abundance ^^Ji 2 is given as usual by the standard calculation of WIMP 
thermal abundance h 2 ). The region where T e > Td is where axinos can dominate the 
universe. This occurs for m„ < 8 TeV in the Tr = 10 10 GeV case. Furthermore, the ratio 
between the T e and To curves gives an approximate measure of the entropy injection from 
axino decay. In frame a), for Tr = 10 6 GeV, axinos never dominate the universe, since 
they decay prior to the point where the axino energy density exceeds that of radiation. 

In frame b). of Fig. |3], we show the same temperatures, but this time for BM2, the 
inoAMSB case with wino dark matter. In this case, the region with mj < 550 GeV is 
excluded by BBN constraints since the axino width is suppressed by the fact that only the 
decay to Z17 is open, where the bino component of Z\ is tiny. When the decay a — > Z27 
turns on around rria ~ 450 GeV, the axino width rapidly increases, and the heavy axino 
scenario becomes BBN-allowed. For Tr = 10 10 GeV, axino domination occurs out to 
rria ~ 10 TeV, while for higher m„ values, the axino decays before neutralino freeze-out. 
For the Tr = 10 6 GeV case, axino domination only occurs in the BBN excluded region. 

In order to see how large the increase in entropy due to axino decay can be, we plot 
in Fig. |] regions of r ranging in value from 1 to > 10 4 in the f a /N vs. Tr plane for 
rria = 1 TeV and a), the BM1 benchmark and b). the BM2 benchmark. The region with 
fa/N > 10 13 GeV is likely BBN excluded since the axino decay temperature To drops 
below ~ 2 MeV. From Fig. |4|, we see that when Tr < T^ cp , and axinos are produced 



via Eq. O, r decreases with increasing f a /N . Also, r increases with increasing Tr due 
to enhanced thermal production of axinos. In contrast, when Tr > T^ cp , the axinos are 
produced in thermal equilibrium and the production rate is independent of f a /N. In this 
case, the r contours increase with increasing f a /N, while they are also independent of Tr. 

6. Axion production from vacuum misalignment 

As already mentioned, the axion will contribute to the total dark matter, so we wish to 
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Figure 3: Plot of To, T e and Tf r versus ma in the BM1 and BM2 benchmark models with 
fa/N = 10 12 GeV. 



review how it is produced. We consider the scenario where the PQ symmetry breaks before 
the end of inflation, so that a nearly uniform value of the axion field 0i = a{x) / (f a /N) 
is expected throughout the universe. As implied by its equation of motion, the axion 
field stays relatively constant until temperatures approach the QCD scale Tqqd ~ 1 GeV. 
At this point, a temperature-dependent axion mass term turns on, and a potential is 
induced for the axion field. At temperature T a the axion field begins to oscillate, filling 
the universe with low energy (cold) axions. The standard axion relic density (via this 
vacuum mis-alignment mechanism) is derived assuming that coherent oscillations begin in 
a radiation-dominated (RD) universe (T a < Tjj or T a > T e ), and its final form is given 
by H H 

fa/N y /6 

10 12 GeV J 



tt s a td h 2 



0.23/(000, 



(6.1) 
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mSUGRA (4525, 275, 0, 1 0, 1 , 1 73.3), m; = 1 TeV 
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Figure 4: Plot of r values in the f a /N vs. Tr plane for a) the FP model and 6) the inoAMSB 
model, with 7715 = 1 TeV. 



where < 9% < vr and f{9i) is the anharmonicity factor 

r / - m 7 /6 

parametrize the latter as f(0i) = In 



Visinelli and Gondolo [34] 
The uncertainty in £l a h 2 from vacuum 



mis-alignment is estimated as plus-or-minus a factor of three. 

However, if the axion oscillation starts during the matter dominated (MD) or the 
decaying particle dominated (DD) phase (Tjj < T a < T e ), the axion relic density will no 
longer be given by Eq. 6T. The appropriate expressions for each of these cases are given 
in Appendix [B[ 



7. Relic abundance of neutralinos 



In Ref. [18], formulae for the relic abundance of neutralinos are derived which include 



the effects of enhancement from axino production and decay and diminution from the 
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re-annihilation effect and entropy dilution. The starting point to evaluate the neutralino 
abundance is the Boltzmann equation 



^ + 3^ i = -H4 i , (7.1) 

where is the neutralino number density, H(t) is the Hubble constant at time t (H 2 = 
p(i)/3M|,), a is the neutralino annihilation cross section, v is the Z\ — Z\ relative velocity 
and (• • •) denotes thermal averaging. At very early times and high temperatures, is 
given by the thermal equilibrium abundance, but as the universe expands and cools, at a 
temperature Tf r the expansion rate H outstrips the annihilation rate n^^(av), and a relic 
abundance of neutralinos freezes out. Thus, we define here the freeze-out temperature by 
the value of T at which 

W% * H{T fr ). (7.2) 
From this condition we can compute the yield variable Y = n^Js: 

Yt{T ir ) = k fr X^L = ffifl. (7.3) 
Z\ s (av)s 

where kf r = 1(3/2) if the freeze-out happens in the radiation (matter) dominated phase. 
If we assume a nearly constant value of (av) (which is appropriate for a higgsino- or wino- 
like Z\ which dominantly annihilates via s-wave reactions), and freeze-out in a radiation- 
dominated universe with H 2 = p ra ^/3Mp and p ra d = 7r 2 g , *T 4 /30, then one obtains 

/r _ (90/yrV(T /r .)) 1/2 

Zi A(av)M P T fr ' y ' 1 

The usual cosmological assumption is that the yield Y~ r is conserved from T = Tf r to 
To, where To is the present day temperature of radiation To = 2.725°K. In this case the 
thermal neutralino relic density is simply 

Zi 45 p c /h 2 Zl z 1 y > 

ft 

However, just as in the axion case, we must re-evaluate Y~ at T = To, and also the freeze- 
out temperature, in the cases of a MD or a DD universe. The various yield expressions are 
contained in Appendix |^. 

7.1 Neutralino re-annihilation at T = Tq 

In the mixed aZ\ DM scenario, neutralinos are produced via axino decay at temperature 
T = Tjj, as well as via thermal freeze-out at T = Tt r . The needed condition for re- 
annihilation at T ~ Tp is that the annihilation rate exceeds the expansion rate at T = Tp: 

-aPW > mM 
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Figure 5: Plot of regions of ma vs. Tr plane where Tf r < T e (right of yellow contour) and where 
no additional neutralino annihilation occurs (brown), in the BM1 (HB/FP) benchmark model. We 
also show regions of entropy generation r. 

where (Tjj) is the total neutralino number density due to thermal (freeze-out) and 
non-thermal (axino decays) production. We can rewrite the above condition as: 

/ - a fr A (90/7rV(T D )) 1/2 
Y~ Zi \t=t d = [Y~ Zi + (T D )) > A{av)MpTD , (7.7) 



where for r > 1 we have Y~ = Ya/r and Y~ r (To) is given by Eq. p.12 . If this condition 



is satisfied, then additional neutralino annihilation will occur at T ~ To- 

The relevant regions are shown in the Tr vs. ma plane for f a /N = 10 12 GeV in Fig. 
|5| for BM1 and in Fig. || for BM2. The magenta, blue and green regions show ranges of 
different entropy dilution factors r, while the brown-shaded region is where there are not 
enough neutralinos produced through axino decays to cause re-annihilation. The region to 
the left of the BBN line is excluded, since To < 2 MeV, and thus we see that for nearly 
all of the allowed parameter space, neutralino re-annihilation effects need to be included 
in our calculations of the neutralino relic density. We also show the (vertical, dotted) line 
where Tjj = Tt r . The region to the right of this line is where Trj > Tt r and the axino 
decay products are expected to thermalize before neutralino freeze-out . 

A simple expression for the diminution of neutralinos from re-annihilation has been 



worked out in Ref. [18] in the sudden-decay approximation. Recasting the Boltzmann 



equation 7.1 in terms of the yield variable 



-^ = - {av) Yls (7.8) 
and integrating from time t = to t gives 

{av)s(T D 



Y,\T) = y-pfc) - M (± - . Y, z \T D) + 



H(Tj 



D 



(7-9) 
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inoAMSB, m 3/2 = 50 TeV, tanp = 10, f a /N = 10 12 GeV 
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Figure 6: Plot of regions of ma vs. Tr plane where Tf r < T e (right of yellow contour) and where 
no additional neutralino annihilation occurs (brown), in the BM2 (inoAMSB) benchmark model. 
We also show regions of entropy generation r. 



where Ys (Tjj) = Y~ r + Y~ , and for Tr > 10 4 GeV, the decay contribution can be huge, 

Z\ Z\ 

so that y~ 1 (To) is small. Then Eq. 7.9 is dominated by the second term and 

Y~ (T) ~ H(T D )/(av)s(T D ) ~ y£(T D ) (7.10) 

Since the latter term is evaluated at T = Tr>, it is much larger than Y~ r (Tf r ), and in fact 
the neutralino abundance turns out to be nearly the freeze-out abundance but evaluated 
at the much lower temperature Tp. 

To exhibit the effect graphically, in Fig. ^ we plot the neutralino yield Y versus Tr for 
a), the BM1 benchmark and b). the BM2 benchmark. In both cases, we take f a /N = 10 12 
GeV and m„ = 1 TeV. The blue curve represents the sum of the thermal neutralino yield 
from freeze-out plus the neutralino yield from axino production and decay, where this term 
dominates for all but the lowest values of Tr in the plot. It increases with Tr as the thermal 
yield of axinos, Eq. 4.3, increases linearly with Tr. For Tr high enough that axinos can 
dominate the universe, the neutralino abundance is diluted by entropy production by the 
same factor as axinos are produced, and thus the curve becomes flat. The green curve is 
the second term of Eq. |7.9| and corresponds to the thermal neutralino yield, Eq. [7.4| , but 
evaluated at Tr> instead of Tf r . The red curve denotes the final neutralino yield, i.e. that 
given by Eq. 7.9. From Fig. f| we see that, for Tr ~ 10 4 GeV, the final neutralino yield 
(Y^ (T)) is close to the naive sum of thermal and non-thermal production (Y^ (Tp)) with 
a small suppression due to the re-annihilation effect. As Tr increases, Y~ increases and 
so does Yg (T), but at a smaller rate, since the re-annihilation becomes more and more 
efficient. Once Tr ~ 10 5 GeV, the annihilation term is efficient enough to re-annihilate 
all neutralinos down to their equilibrium value (Y"~ r (Trj)). Thus, for Tr > 10 5 GeV, the 

neutralinos injected from axino decays can effectively thermalize to Y~ r (Ti)) and the final 

Zi 
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mSUGRA (4525, 275, 0, 10, 1, 173.3), f,/N = 10 12 GeV, m;= 1 TeV 

10~ 4 
10" 5 

itr 6 

10' 7 

10" 8 

10" 9 
10" 10 
10"" 

10 4 10 5 10 6 10 7 10 8 10 9 10 10 10 11 10 12 

T R [GeV] 

inoAMSB, m 3/2 = 50 TeV, tanp = 1 0, f a /N = 1 12 GeV, m~ = 1 TeV 
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Figure 7: Plot of yield Y versus T R for f a /N = 10 12 GeV and m a = 1 TeV for a) the FP model 
and b) the inoAMSB model. 

neutralino yield becomes independent of the initial axino abundance. Since Tjj is usually 
far below Tf r , Eq. [7.4| tells us that the neutralino yield Y^(T) ~ Y~ (To) will likely far 

exceed the thermal yield evaluated at Tf r , Y~ r (Tf r ). 

For benchmarks BM1 and BM2, the region to the right of the yellow contour and to 
the left of the vertical blue line denotes where, in Figs. || and ^, we have To < Tf r < T e . In 
this case, the thermal neutralino abundance must be calculated in a MD or DD universe. 
To illustrate, we show the various neutralino yields in Fig. |§] for benchmark BM1 with 
rria = 5 TeV. The pink curve denotes the thermal neutralino abundance, which is constant 
versus Tr until r > 1, where the abundance becomes diluted by 1/r and begins to decrease. 
At Tr ~ 2 x 10 10 GeV, T e begins to exceed Tf r , the neutralino freeze-out happens in a MD 
universe, and suffers an increase. The corresponding freeze-out temperature is denoted 
by the azure colored curve, and the right-hand y-axis: we see that T* r drops from ~ 3.8 
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Figure 8: Plot of yield Y versus T R for f a /N = 10 12 GeV and m~ a = 5 TeV for the BM1 (HB/FP) 
model. The dashed gray line shows where r = 1 and the yellow line shows where Tr = T c i cp - 

GeV to ~ 3.5 GeV as we move from a RD to a MD universe. From Fig. |8] we see that 
the change in Y~ r is inconsequential to the final neutralino yield (Y^(T)), since this is 
determined by the re-annihilation term (green curve), and is insensitive to the thermal 
neutralino abundance, at least in this case. 

8. Algorithm for DM abundance in aZ\ DM scenario 

In this section, we list our algorithm for evaluating the relic density of mixed aZ% CDM in 
the PQMSSM. We proceed as follows. First, extract an effective value of (av) by matching 
Eq. |^4| onto the fl^ h 2 result from the IsaReD subroutine of Isajet (the effective (av) 
value is shown in Table 111). Then: 



If To <Tbbn (where Tbbn is taken as 2 MeV), regard as BBN excluded) 35] 
If T D > T fr , then n^h 2 = fi| d /i 2 , given by the IsaReD output. 
If To < Tf r and Tp > Tbbn, then 

— If r > 1: (Case of axino domination with T e > Tp) 



* If Z\ re-annihilate (Eq. |7\7] satisfied), then calculate using Eq. /\£ with 



Y~ (r D ) = Yl T + YS , where Yl r is given by Eq. C13 and YS = Y~Jr. 



* If Z\ does not re-annihilate (Eq. |7T^ not satisfied), then Y~ = Y~ r + Y~ , 



where Y~ is given by Eq. |C.12| and Y~ =Y%/r. 

Z\ Z\ 



If r < 1: (axino non-domination) 
* If neutralinos re-annihilate (Eq. |7.7| holds), then yield given by Eq. 7J} with 



Yl r given by Eq. \l.4 and Y% =Y~ a . 
Z\ — Z\ 
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If neutralinos do not re-annihilate (Eq. |7.7| does not hold), then yield given 
by Yl r + Kj with Yl r given by Eq. O and Y% =Y- a . 

Z\ Z\ Z\ 1 1 Z\ 



• Now add in axion contribution to relic abundance Eq. B.14 

• Final dark matter abundance is given by: Q a ^^h 2 = Slj^/i 2 + Q a h 2 . 

9. Numerical results 

In this section, we present results for the relic abundance of mixed aZ\ CDM for benchmarks 
BM1 and BM2 in the PQMSSM. Our first results are shown in Fig. ||, where we plot the 
neutralino abundance, il^h 2 , the axion abundance, O a /i 2 , and their combination, Q^^h 2 , 
versus m„ for f a /N = 10 12 GeV, with Tr = 10 10 GeV. We take the initial axion field value 
Oi = 0.498 for the BM1 benchmark case shown in frame a), and we take 0i = 0.675 for the 
BM2 benchmark in frame b). These values tune the total dark matter abundance to the 
WMAP value for the case where Td > TV r . From frame a)., we find the region to the left 
of the dashed gray line is excluded by BBN constraints on late decaying neutral particles, 
since Td < 2 MeV. For m a values just beyond the BBN constraint, the naive expectation 
is that the neutralino abundance is determined by the thermal axino production rate, and 
indeed Fig. suggests the abundance Q% h 2 ~ 10 7 . Instead, the actual abundance is 
several orders of magnitude below this, but still far above the measured DM value. In 
this case, the large thermal axino production rate is followed by decays to neutralinos at 
T = Tjj. As relic neutralinos fill the universe they proceed to annihilate, so that their 



final abundance is determined by Eq. 7.9. Since in this region ft^/i 2 ~ l/T D , and since 



3/2 , 

Td ~ rn a , we find the neutralino abundance decreasing with increasing m a . The kink 
in frame a), at m„ ~ 900 GeV occurs due to turn-on of the a — )• gg decay mode, which 
increases Tq, thus decreasing ^^/i 2 even further. The jog at mj ~ 1.6 TeV is caused by 
a change in the g* value due to the addition of quark and gluon degrees of freedom after 
the QCD phase transition. While ^^ ± h 2 is decreasing with increasing m^, it reaches 0.11 
at rria ~ 6 TeV and continues dropping until To exceeds Ty r . At this point, the thermal 
Z\ abundance assumes its traditional value of £l s ~ d h 2 ~ 0.05 since now axinos decay before 
freeze-out. For rria ^ 8.5 TeV, the CDM is a nearly equal mix of axions and neutralinos. 
While Tq > 1 GeV (for high mj), the axion abundance assumes the form as given in Eq. 



6.1. However, if considering lower values of m a , then r > 1 and Td < 1 GeV, so that 
the axion abundance is diluted by entropy production from axino decay. As m a increases, 
the axion dilution becomes lessened until Td exceeds ~ 1 GeV, after which the axion 
abundance remains fixed. 

In frame b)., for the inoAMSB model, we again see that the neutralino abundance is 
large at small m a : ^l^Ji 2 ~ 10 for m a ~ 600 GeV. It decreases steadily with m a due to the 
increasing Td- The curve does reach a point where the neutralino abundance matches the 
WMAP7 measured value, at m a ~ 1.5 TeV. In fact, the neutralino production via axino 
decay and subsequent re-annihilation effect offers an elegant means to enhance the WIMP 
relic density in SUSY models with higgsino or wino-like WIMPs. This mechanism offers an 
alternative [26] to the neutralino abundance enhancement via moduli decays as proposed 
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inoAMSB, m 3/2 = 50 TeV, tanp = 1 0, f a /N = 1 1 2 GeV, T R = 1 1 GeV, 8i = 0.675 
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Figure 9: Plot of neutralino and axion relic densities Vth 2 versus m~ a for f a /N — 10 12 GeV and 
T R = 10 10 GeV for a) the HB/FP model and b) the inoAMSB model. 



by Morio and Randall |3q] for AMSB, and as proposed by Kane et al. to explain the 
Pamela/Fermi cosmic ray anomalies [37]. By the time Tp exceeds Tf r , the axion abundance 
has assumed the value given by Eq. and % ~ 1CT 3 , as given m Table @. 



In Fig. 10, we show the mixed aZ\ abundance versus f a /N for fixed m„ = 1 TeV 
and T R = 10 10 GeV for a). BM1 and b). the BM2 benchmarks. In the HB/FP case 
of frame a)., we see that for low f a /N, the axino width Ta is large, and Tjj exceeds 
Tf r , so that axinos decay before freeze-out and the neutralino relic density assumes its 
standard value of Q, s ~ d h? = 0.05 in this case. Meanwhile, the axion density is extremely 
small due to the low value of f a /N . As f a /N increases, the axion abundance naturally 
increases, while the neutralino abundance remains constant until around f a /N ~ 10 10 
GeV, where To falls below Tf r . For higher f a /N values, To continues to fall and since 
■-' -■ Tp , the neutralino abundance steadily increases. At the point where r = 1 is 
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Figure 10: Plot of neutralino and axion relic densities fth 2 versus f a /N for rria = 1 TeV and 
T R = 10 10 GeV for a) the HB/FP model and b) the inoAMSB model. 



reached (f a /N ~ 2 x 10 10 GeV), entropy injection from axino decay causes a small decline 
in the otherwise steadily increasing axion abundance. The general behavior in frame b). 
for the inoAMSB model is similar to that of frame a). Nowhere in these two plots does 
the axion abundance exceed the neutralino abundance. This is merely a reflection of our 
choice of m„ = 1 TeV and Of, for higher values of m„, the value of Tjj increases, resulting 
in a decrease of £1^ h? when re-annihilation dominates the neutralino relic abundance. 

In Fig. 11, we show the mixed dark matter relic abundance versus Tr for fixed f a /N = 
10 12 GeV and fixed rria = 1 TeV, for benchmarks BM1 and BM2. In this case, Tp is fixed 
throughout the plots, and so h 2 is nearly constant everywhere except at low Tr ~ 
10 4 GeV, where thermal axino production is somewhat suppressed, and fewer neutralinos 
are produced at To to enter the re-annihilation process. Since f a /N is fixed, the axion 
abundance is also constant throughout much of the plot. At Tr > 10 8 — 10 9 GeV, we 
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Figure 11: Plot of neutralino and axion relic densities Qh? versus Tr for rria = 1 TeV and 
f a /N = 10 12 GeV for a) the HB/FP model and b) the inoAMSB model. 



enter the region where axinos can dominate the universe (r > 1), and entropy production 
from axino decay diminishes the axion abundance. For even higher values of > ^dcpi 

the 

axino production rate becomes independent of Tr, and the entropy injection ratio r becomes 
constant with Tr. While the neutralino abundance dominates the axion abundance in these 
frames, again, this is just a reflection of the value of mj chosen; for higher m„, To will 
increase, leading to a diminution of O^/i 2 . 

We have seen that over most of parameter space with Tp < Tj r , Qg^ 2 ~ l/Tp ~ 

^ i<2 3/2 

T- ~ (f a /N)/m~ , with little dependence on Tr. Hence, a good way to display the relic 
density of dark matter in the mixed aZ\ CDM scenario is to display it in the vs. (f a /N) 
plane for benchmarks BM1 and BM2. This plane is shown in Fig. |l^a). for BM1 and frame 
b). for BM2. Here, we take 6i = 0.498 so as to normalize the relic density ^ a ^^ 2 to the 
measured value 0.1123 when Tjj > Tj r and f a /N = 10 12 GeV. The black contour denotes 



the line where Td = Tf r : below and right of this contour, the neutralino relic density 
is given by its usual thermal abundance, which is ^ Zl h 2 = 0.05 for the BM1 case in 
frame a). In this region, the axion abundance increases with increasing / Q /iV, so that 
fi^/i 2 = 0.1123 at f a /N = 10 12 GeV by design, with a roughly even admixture of mixed 
higgsino and axion dark matter in the narrow azure-shaded region. In the region to the left 
of the Tq = Tf r contour, the neutralino abundance rapidly increases, and we have regions 
of dominantly WIMP CDM. 

In frame b). for BM2, we again adjust 9i so that the total relic density equals the 
measured value for Td > Tf r with f a /N = 10 12 GeV. We see qualitatively similar behavior 
as in frame a). The azure region for Td > Tf r (right-side of plot) has fi^r ~ 0.02, and so 
is axion-dominated, while the azure region for Td < Tf r is wino-dominated. In the region 
with Td < Tf r , the relic abundance of winos rapidly increases as we move to smaller 
or larger f a /N values. 

From the above results for the benchmarks BM1 and BM2, we see that the aZ\ scenario 
can be classified into two main cases: 



• A): decoupled axino (Td > Tf r ) 

• B): axino enhanced DM (Td < T* r ) . 

Case A) happens for high Td values, which are obtained at low f a and/or high ma, as 



seen in Figs. |9j and 10. In this scenario, the axino has no effect on the DM relic density, 
which can be a mixture of axions and neutralinos. Since the axion mis-alignment angle 
(6i) can always be adjusted so that Q a h 2 = 0.1123, there is no lower bound on VL^J% 2 . 
Nonetheless the neutralino relic density must still satisfy: 

n^h 2 = Sl s gh 2 < 0.1123 (T D >T fr ) (9.1) 

where Q~ d h 2 is the standard neutralino freeze-out relic density in the MSSM, since there 
is no axino dilution or contribution in this case. Therefore, in Case A, any MSSM model 

is composed of axions. 

For Case B), Td < Tf r , which is obtained at high f a and/or low mg. As shown in 
Figs. ^ and [||, for most of the parameter space considered here, the neutralino relic density 



satisfying Eq. |9.l| is allowed. For models where Q s ~ h < 0.1123, the remaining of the DM 



is dominated by the annihilation term in Eq. |7.9| . In this case the relic density can be 
approximated by: 

% h 2 ~ nfh 2 x (9.2) 
Assuming Tf r ~ m^/20 and using Eqs. 5.2 and 3.2, we obtain: 

n~ h 2 ~ 25 x n-^h 2 ( m ^ ) ( fa/N ) f 1 -^] 3/2 (i - ^l) ~ 3/2 (9 3) 

U z h -2bxU Zi h ^ 100 GeV A 10 12 GeV J ^ m~ a ) V m 2 J ^ 
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Figure 12: Regions of neutralino plus axion relic density Sl^/i 2 in the mg vs. f a /N plane for 
T R = 10 10 GeV for a) the FP model and b) the inoAMSB model. The unshaded (white) regions 
are excluded by BBN bounds since Tr < 2 MeV. 



where we assumed m„ > nig. Now imposing the DM relic density constraint (Eq. |9.l|) , we 
obtain: 

w < 4 „ ur. (~) (• - 3) 3/2 < r » < ^ 

(9.4) 

Therefore, in this case, the MSSM relic density has to be considerably suppressed in order 
to satisfy the above bound. Although the bound decreases with f a and increases with ma, 
for sufficiently low f a and/or high mg, then Tjj > Tf r and the bound in Eq. |9.l| must 
be used instead. In the case where Tp < Tf r , the DM will likely be composed mainly of 



relic neutralinos, unless Q~ d h 2 is much smaller than Eq. 9.4. We also point out that the 



approximate bound in Eq. 9.4 is a conservative one, since, for < mn, the bound would 
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Figure 13: Neutralino, axion and total relic density for the inoAMSB model versus m 3 / 2 with 
tan/3 = 10 and \i > and for nia — 2 TeV and 8i = 0.675. We also show the value of r by the 
yellow curve, and the right-side axis. 



be more strict. 



To see how these conclusions depend on the SUSY spectrum, we show in Fig. 13 the 
neutralino, axion and summed relic abundance for the inoAMSB model versus TO3/2 for 
tan/3 = 10 and \x > 0. The results hardly change with varying tan/3 or fi. We also take 
f a /N = 10 12 GeV, m- a = 2 TeV, 9 t = 0.675 and T R = 10 10 GeV. As m 3/2 increases, all the 
sparticle masses increase as well and so T^ decreases. Once m 3 / 2 ~ 50 TeV, rrig ~ m„/2 
and the bound in Eq. 9.4 can no longer be satisfied. The ratio of entropy injection r 
is shown as the yellow curve, against the right-hand y-axis. Since r = AmgYa/'STD, and 
Td ~ r^ 2 , Td decreases with increasing m 3/ / 2 , and the entropy ratio increases. This leads 
to a dilution of the axion relic density as shown in the plot. The jog in the curves around 
m 3/2 ~ 55 TeV occurs due to a change in the degrees of freedom g*. 



10. Effect of saxion production and decay on relic abundance 

We have seen so far that the relic neutralino abundance may be enhanced beyond usual 
expectations in the mixed aZ\ scenario if Td < Tf r . However, we have so far neglected a 
mandatory element of the axion supermultiplet: the spin-0, i?-parity even saxion field sp8|. 
Saxions may be produced thermally in the early universe, either in thermal equilibrium 
for Tr > Tdcp or via radiation and decay for Tr < Td cp - Saxions can also be produced via 
coherent oscillations. 

In an analagous manner to the axino case, the saxion field can dominate the energy 
density of the universe at early times, if its decay temperature (Td s ) is smaller than the 
temperature (T 6s ) at which its energy density overcomes the radiation and axino energy 
densities. The saxion can decay into gluons and gluinos (and perhaps axions), depending 
on its mass and its (model dependent) couplings. As in the axino case, the saxion decay 
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will inject entropy at T = Td s , diluting the axino (if Tjj < Td s ), neutralino (if Tf r > Tn s ) 
and axion (if T a > Td s ) relic densities. Furthermore, if the s — >■ gg branching ratio is 
considerable, saxion decays will also inject neutralinos through gluino cascade decays. In 
the case where a high rate of neutralino injection occurs after freeze-out, one must again 
consider a second possibility of neutralino re-annihilation which may enhance the neutralino 
abundance. 

Since the saxion lifetime is comparable to the axino lifetime, it is possible that both 
saxions and axinos may co-dominate the universe. This makes a simplistic analysis of 
neutralino abundance difficult using the approach of this paper. To make matters worse, 
gravitinos may be produced thermally at significant rates for high enough Tr, and enjoy 
decays at time scales comparable to heavy axinos and saxions. Thus, gravitino production 
and decay may also enhance or diminish the neutralino abundance. 

The proper treatment of such intricately coupled effects is best made by numerical 
solution of the coupled Boltzmann equations. This type of treatment has been initiated in 



Ref. |3J|, and will be reported at a future date. 



11. Conclusions 

In this paper, we have presented results of calculations of the dark matter abundance in 
supersymmetric models wherein the strong CP problem is solved by the Peccei-Quinn 
mechanism, and in which the neutralino is the LSP, so that dark matter consists of an 
axion/neutralino admixture. Since the aZ\ CDM scenario boosts the dark matter abun- 
dance beyond the usual thermal neutralino production rates, we have presented results for 
two models that typically yield an underabundance of thermal neutralino dark matter: the 
HB/FP region of mSUGRA with a mixed higgsino-like neutralino and AMSB-type models 
with a wino-like neutralino. Our final results depend mainly on the temperature Tr> at 
which heavy axinos finish their cascade decays to neutralinos. 

In the case where Tp > Tf r , the neutralino abundance is given by its usual thermal 
abundance. In the case where there is an underabundance of neutralinos compared to 
the measured dark matter abundance, then the remainder can be comprised of axions. In 
the case of our two benchmark models, in the HB/FP region, we would obtain a nearly 
equal mixture of axions and mixed higgsino-like neutralinos, while in the BM2 case of the 
inoAMSB model, we would obtain a case with mainly axion CDM, along with a small 
admixture of wino-like neutralinos. The case of BM1 provides an instance where both 



WIMP and axion signals [43] could ultimately be found at dark matter detectors. In the 
case of BM2, only an axion signal might be discovered. The direct detection rates for 
wino-like WIMPs has been presented in Ref. | ]26| |. These projected rates would have to be 
scaled down by a factor of ~ 70 since in this case wino-like WIMPs would only comprise 
~ 1/70 of the total dark matter abundance. 

In the case where Td < Tf r , heavy axino decays in the early universe will inject 
additional neutralinos and entropy into the cosmic soup after neutralino freeze-out. The 
main effect here is that neutralino re-annihilation takes place at T = Tp, and one obtains 
a neutralino abundance as if the freeze-out temperature were replaced by Tjy. In this case, 
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^Zi^ 2 ~ V^£> instead of l/7/ r , and since Tp> < Tf r , one obtains a greatly enhanced 
abundance of neutralinos beyond the usual expectation. The neutralino "production by 
decay and re-annihilation" mechanism can thus lead to enhanced production of wino-like or 
higggsino-like neutralinos- which naively give rise to an underabundance of dark matter- 
so that Zis might comprise nearly all the abundance of dark matter. Also in this case, if 
T e > Td, then significant entropy production by decaying axinos can diminish the axion 
abundance. In this case, one might expect mainly wino-like or higgsino-like dark matter, 
with a small admixture of axions. Thus, the enhanced neutralino production via axino 
decay mechanism offers an alternative means to allow wino- or higgsino-like neutralinos 
to comprise the bulk of dark matter. In many respects, this mechanism may be preferred 
over the possibility of multi-TeV scale moduli field decay []36|, [3^1 , since it also allows for a 
solution to the strong CP problem. 
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A. Expansion rate of early universe 



Here we briefly review the cosmology of an early axino dominated universe and in the 
next Section present the expressions for the axion relic density and neutralino yield used 
in Sees. 

First, we define several temperatures: 

• T e : temperature when the universe becomes axino (matter) dominated, 

• Ts- temperature at which entropy injection due to axino decay starts, 

• Tp>: temperature at which entropy injection due to axino decay ends. 

A.l Matter dominated phase: T s <T <T e 

In this phase, the universe is axino dominated, for which: 

R? R? 
Pa = Po-^fexp(-(i - t e )T- a ) « po-^f (A.l) 

H(T) = , [Jj£ ~ . ( ^) ^ (A.2) 

V ' V 3MJ, V 3Mp \ RJ V ; 

where R e = R(T e ) and r„, pg, Y„ and rria are the axino width, energy density, yield and 
mass, respectively. Here, Mp = Mpi/-v8n (i.e. Mp is the reduced Planck mass), and we 
will use g*s = 9* since the temperatures we consider are all above Tbbn — 2 MeV. By 
definition we have: 

7T 2 4 4 

PO = ^9*( T e) T e and T e = ^0^0,- (A. 3) 
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B. Axion oscillation 



The axion field starts to oscillate when 2 

3H(T a ) = m a (T a ) 
where the temperature-dependent axion mass is given by 

m a , if T < A 



m a (T) 



m a b(±) 4 \ if T > A 



(B.l) 



(B.2) 



with b = 0.018, A = 0.2 GeV and m a = 6.2 x 10 3 // Q . Due to its temperature dependent 
mass, after oscillation begins, the axion energy density obeys 



Pa{T)R{Tf 



constant. 



(B.3) 



m a (T) 

If T a < Td, then the axion field starts to oscillate after the axino has decayed. In this 



case, the axion relic density is given by the standard expression [34, 42]: 



m td h 2 



/ , \ 3/2 

9.23 X 10- 3 ^/(^)^yT74 (i^J , if fa > fa 

7/6 



1.32 0ff(6i 



fa 



(B.4) 



, if fa < fa 



where f a = 9.9 x 10 16 GeV, f(6i) = ln( „|, 2 ) 7 ^ 6 and the standard oscillation temperature 

i I 

is given by 

1.23xltf« ;: ^ FI (l^) 1/a > if/ a >/ a 



8.71 x 10 



-l l 



T) ' ,if/.</. 



1/6 



(B.5) 



If instead T e < T a , then the axion density is diluted by the entropy ratio r so that 



Slab? 



1 



x n s „ td h 2 



(B.6) 



If Td <T a < T e , the axion can start to oscillate in the matter dominated phase (MD) 
or the decaying dominate phase (DD). The relic densities for each case are: 



• Matter dominated (T$ < T a < T e ): 



7.5 x 10- 5 e1f{Bi)T D I -1* 



10 1 



, if fa > f, 



MD 



1 



1-4 2 f(9i) gt{Ta)4/11 y w 



14/11 



r, 



l d ;f f ^ fJS 
47TTj 11 J a ^ J a 



MD 



with 



and T a given by 



fa ID = 7.6 x 10 



17 



Vd*{Ta)T e 



(B.7) 



(Bi 



2 Here, we follow much of the notation given by Visinelli and Gondolop^]. 
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-iMD 



6.1 x 10 2 



i 10 1 



8.6 x 10 



-i 



y/9*(T a )T e U 



1 10 1 



2/3 



, ^ fa > f 



AID 



2/11 



y fl .(T a )T e /a 

Decaying particle dominated (Td <T a < Tg)p^]: 



if/«</a J 



7.5 x 10- 5 ^ 2 /(^)T D (^J , if f a > f° D 
1 72 P 2 ff6M g * (7b)1/4 ^ 2 fJ^ 3/2 if f / f£"- ; 



ZShfr) ,iffa<f a U 



with 



and T a given by: 



fDD = 5 6g x 1Q 



20 



Vg*(Tp) 

9*{Ta) 



T, 



D 



^DD 



0,1 xlO 2 ^^ 2 ) ,if/ a >/r 

1/8 



9.0 x 10 



-J / y/g*{Tp) 10 12 rp2 \ if f , iDD 

9*{T a ) fa 1 D) > 11 /a < Ja 



Matching both solutions we have: 



fa A ^ <Q MTa) V22 1 



12 



10 

fa 

To 12 



5/22 



> 0.8 



gjjaf^ i 
9.(T D y/* TDJ t/u 



MD case 



DD case 



or in terms of T a : 



T a > T s -> MD case 
T a < Ts — > DD case 



(B.9) 



(B.10) 



(B.ll) 



(B.12) 



where 



To summarize: 



9*(T D ) 



1/5 



,3* (7a 



rr-i r7-:4 



n n h 2 



U a td h 2 /r , if T e < T a 
fif D /i 2 , if T 5 < T a < T e 
^/j 2 , if Td < T a < T 5 
[ n s a td h 2 , if T a < Td 

where r is the entropy injection ratio as usual. 



(B.13) 



(B.14) 
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C. Neutralino yield 



The neutralino will decouple from the thermal bath when 

(av)n 2i (T fr ) = H(T fr ) (C.l) 

where 



/mT\ 3/2 



-m/T 



(C.2) 



In a radiation dominated universe, the neutralino yield is given by: 



while in a matter dominated universe: 

As in the axion case, the neutralino can freeze-out before the universe becomes matter 
dominated (Tj r > T e ), during the MD phase (T5 < Tf r < T e ), during the DD phase 
(To < Tf r < Ts) or during the radiation dominated phase (Tf r < Tjj). The neutralino 
yields for each of these scenarios are listed below. 



Standard case (Tf r < Tp): 



r tm _ W^f'f 2 ,c,) 



Zi v J " A{av)M P Tf 
where the freeze-out temperature is given by 

3 v / 5(^> M P m% 2 



Tf f d = m g /ln[ Zl ]. (C.6) 



MD case (T s < T fr < T e ): 



V MD(rp \ _ ^ sfrj^pMZ^ Tp 

y~ (T D ) - -y- (T /r (C7) 



where the freeze-out temperature T^ D is given by 



3V5(av)M P mf 
T¥ r D = m~/\n[ — %-\. (C.8) 



DD case (To <T fr .<T s ): 
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YP D (T D ) = ^(T^ 



ddn 91 /2 (Td) ( T D 



2-* ' g y\Tf°) \Tf° 
where the freeze-out temperature is given by 



(C.9) 



T 



DD 



3V5(<Jv)Mpm 3 J 2 gl /2 (T D )Tl 



Case when Tf r > T e : 



(CIO) 



Y Si (T D ) = Yj£(Tp)/r = Y§P(Tp) x g. 



std I rpstd\ 



To summarize: 



D 



' Y^(Tf r d )/r , if T e < T fT 
YfD , if T 5 < T fr < T e 
YP D , if T D < Tf r < T s 
Y f(Tf r d ) , if Tf r < T D 



(C.ll) 



(C.12) 



where 



g*{Tp) 



1/5 



-TfjT e 



(C.13) 



References 

[1] H. Baer and X. Tata, Weak Scale Super symmetry: From Superfields to Scattering Events, 
(Cambridge University Press, 2006). 



[2 
[3 
[4 

[5 
[6 
[7 
[8 
[9 

[io; 
[ii 



[12 



H. Baer, V. Barger, A. Lessa and X. Tata, \J. High Energy Phys. 1006 (2010) 102 
G. 't Hooft, \Phys. Rev. Lett. 37 (1976jJ . 



For a recent review, see P. Sikivie, hep-ph/0509198; M. Turner, Phys. Rept. 197 (1990) 67; 



J. E. Kim and G. Carosi, Rev. Mod. Phys. 82 (2010) 557 



For a review, see J. E. Kim, Phys. Rept. 150 (1987) 1 



For a comprehensive discussion, see M. Pospelov and A. Ritz, Annals Phys. 318 (2005) 119. 



R. Peccei and H. Quinn, \Phys. Rev. Lett. 38 (1977) 1440| and \Phys. Rev. D 16 (1977) 1791 . 
S. Weinberg, \Phys. Rev. Lett. 40 (1978) 223| ; F. Wilczek, \Phys. Rev. Lett. 40 (1978) 279 . 
J. E. Kim, \Phys. Rev. Lett. 43 (1979) 103fc M. A. Shifman, A. Vainstein and V. I. Zakharov, 



NucL Phys. B 166 (1980) 493 



M. Dine, W. Fischler and M. Srednicki, \Phys. Lett. B 104 (1981) 199J A. P. Zhitnitskii, Sov. 
J. Nucl. 31 (1980) 260. 

Phys. Rev. D 18 (1978) 1829] and \Phys. Rev. 



D. Dicus, E. Kolb, V. Teplitz and R. Wagoner, 



D 22 (1980) 839| ; for a review, see G. Raffcldt, |hep-pti/06 11350 



H. P. Nilles and S. Raby, Nucl. Phys. B 198 (1982) 102; J. E. Kim and H. P. Nilles, Phys 



Lett. B 138 (1984) 15C. J. E. Kim, Phys. Lett. B 136 (1984) 378 



- 29 - 



[13 

[14 
[15 

[16 

[17- 
[18 
[19 
[20 
[21 



P. Svrcek and E. Witten, \j. High Energy Phys. 0606 (2006) 051 



K. Rajagopal, M. Turner and F. Wilczck, Nucl. Phys. B 358 (1991) 447 



L. Covi, J. E. Kim and L. Ros zkowski, \Phys. Rev. Lett. 82 (1999) 4180| ; L. Covi, H. B. Kim 



J. E. Kim and L. Roszkowski, [J. High Energy Phys. 0105 (2001) 033 



For recent reviews of axino dark matter, see F. Steffen, Eur. Phys. J. C 59 (2009) 557; L. 
Covi and J. E. Kim, iVeto J. Phys. 11 (2009) 105003. 



H. Baer, A. Box and H. Summy, \J. High Energy Phys. 0908 (2009) 080 



K-Y. Choi, J. E. Kim, H. M. Lee and O. Seto, Phys. Rev. D 77 (2008) 123501 



E. Komatsu et al. (WMAP collaboration), arXiv:1001.4538 (2010). 

For a recent review, see R. Arnowitt and P. Nath, arXiv:0912.2273 (2009) 

K. L. Chan, U. Chattopadhyay and P. Nath, \Phys. Rev. D 58 (1998) 096004 J- Feng, 



K. Matchev and T. Moroi, Phys. Rev. Lett. 84 (2000) 2322 and Phys. Rev. D 61 (2000) 



075005 ; see also H. Baer, C. H. Chen, F. Paige and X. Tata, \Phys. Rev. D 52 (1995) 2746j 
and \Phys. Rev. D 53 (1996) 624l|; H. Baer, C. H. Chen, M. Drees, F. Paige and X. Tata, 



Phys. Rev. D 59 (1999) 055014; for a model-independent approach, see H. Baer, 



T. Krupovnickas, S. Profumo and P. Ullio, J. High Energy Phys. 0510 (2005) 02C 



22] L. Randall and R. Sundrum, Nucl. Phys. B 557 (1999) 79; G. Giudice, M. Luty, H. 



Murayama and R. Rattazzi, J. High Energy Phys. 9812 (1998) 027 



[23 
[24" 



[25 



[26 

[27; 

[28 

[29 
[30 
[31 
[32 
[33 



[34; 
[35 



S. P. de Alwis, \j. High Energy Phys. 1003 (2010) 078] . 



H. Baer, S. P. de Alwis, K. Givens, S. Rajagopalan and H. Summy, J. High Energy Phys\ 
1005 (2010) 069|; H. Baer, S. P. de A lwis, K. Givens, S. Rajagopalan and W. Sreethawong, [7] 
High Energy Phys. 1101 (2011) 005] , 



H. Baer and A. Box, Eur. Phys. J. C 68 (2010) 523; H. Baer, A. Box and H. Summy, J. High 



Energy Phys. 1010 (2010) 023 



H. Baer, R. Dermisek, S. Rajagopalan and H. Summy, JCAP1007 (2010) 014. 



ISAJET, by H. Baer, F. Paige, S. Protopopescu and X. Tata, tiep-ph/0312045; see also 



H. Baer, J. Ferrandis, S. Kraml and W. Porod, \Phys. Rev. D 73 (2006) 01501C| . 
F. Paige, S. Protopopescu, H. Baer and X. Tata, |hep-pri/0312045| ; 



http://www.nhn.ou.edu/~isajet/ 



H. Baer, C. Balazs and A.Belyaev, J. High Energy Phys. 0203 (2002) 042 



A. Brandenburg and F. Steffen, JCAP0408 (2004) 008. 



A. Strumia, J. High Energy Phys. 1006 (2010) 036 



R. J. Scherrer and M. S. Turner, \Phys. Rev. D 31 (1985) 681 



L. F. Abbott and P. Sikivie, Phys. Lett. B 120 (1983) 133; J. Preskill, M. Wise and F 



Wilczck, Phys. Lett. B 120 (1983) 127; M. Dine and W. Fischler, Phys. Lett. B 120 (1983) 



137; M. Turner, Phys. Rev. D 33 (1986) 889 



L. Visinelli and P. Gondolo, Phys. Rev. D 80 (2009) 035024 



K. Jedamzik, Phys. Rev. D 70 (2004) 063524 and Phys. Rev. D 74 (2006) 10350£ 



- 30 - 



[36] T. Moroi and L. Randall, \Nucl. Phys. B 570 (2000) 455fc B. Acharya, G. Kane, S. Watson 
and P. Kumar, \Phys. Rev. D 80 (2009) 083529 . 

[37] P. Grajek, G. Kane, D. Phalen, A. Pierce and S. Watson, \Phys. Rev. D 79 (2009) 043506^ G. 



Kane, R. Lu and S. Watson, Phys. Lett. B 681 (2009) 151| . D. Feldman, G. Kane, R. Lu and 
B. Nelson, \Phys. Lett. B 687 (2010) 364 



38] J. E. Kim, Phys. Lett. B 67 (1991) 3465; D. H. Lyth, Phys. Rev. D 48 (1993) 4523; S 



Chang and H. B. Kim, [Phys. Rev. Lett. 77 (1996) 591 ; M. Hashimoto, K. Izawa, M. 
Yamaguchi and T. Yanagida, Phys. Lett. B 437 (1998) 44j; T. Asaka and M. Yamaguchi, 



115029 



Phys. Rev. D 59 (1999) 125003t J. Hasenkamp and J. Kersten, phys. Rev. D 82 (2010) 



[39] 
[40] 
[41] 



H. Baer, S. Kraml, A. Lessa and S. Sckmen, arXiv:1012.3769. 
H. Baer, A. Lessa, S. Rajagopalan and W. Sreethawong, to appear. 
J. McDonald, 



Phys. Rev. D 43 (1991) 1063; C. Pallis, Astropart. Phys. 21 (2004) 689 



G. Lazarides, C. Panagiotakapolous and Q. Shafi, Phys. Lett. B 192 (1987) 323;G. Lazarides, 
R. Schacfcr, D. Seckel and Q. Shafi, \Nucl. Phys. B 346 (1990) 193j. 



[42] L. Visinelli and P. Gondolo, [Phys. Rev. D 81 (2010) 063508 

[43] L. Duffy et al.,\Phys. Rev. Lett. 95 (2005) 091304] and \Phys. Rev. D 74 (2006) 012006 



for a 



review, see S. Asztalos, L. Rosenberg, K. van Bibber, P. Sikivie and K. Zioutas, Ann. Rev 
Nucl. Part. Sci. 56 (2006) 293| . 



- 31 - 



